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Abstract. We study the pseudospectrum of a class of non-self adjoint differential 
operators. Our work consists in a detailed study of the microlocal properties, which 
rule the spectral stability or instability phenomena appearing under small pertur- 
bations for elliptic quadratic differential operators. The class of elliptic quadratic 
differential operators stands for the class of operators defined in the Weyl quantiza- 
tion by complex-valued elliptic quadratic symbols. We establish in this paper a simple 
necessary and sufficient condition on the Weyl symbol of these operators, which en- 
sures the stability of their spectra. When this condition is violated, we prove that 
it occurs some strong spectral instabilities for the high energies of these operators, 
in some regions which can be far away from their spectra. We give a precise geo- 
metrical description of them, which explains the results obtained for these operators 
in some numerical simulations giving the computation of "false eigenvalues" far from 
their spectra by algorithms for eigenvalues computing. 

Key words. Spectral instability, pseudospectrum, semiclassical quasimodes, non- 
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1.1. Miscellaneous facts about pseudospectrum. In recent years, there has been 
a lot of interest in studying the pseudospectrum of non-selfadjoint operators. The 
study of this notion has been initiated by noticing that for certain problems of sci- 
ence and engineering involving non-selfadjoint operators, the predictions suggested by 
spectral analysis do not match with the numerical simulations. This fact lets thinking 
that in some cases the only knowledge of the spectrum of an operator is not enough to 
understand sufficiently its action. To supplement this lack of information contained 
in the spectrum, some new subsets of the complex plane called pseudospectra have 
been defined. The main idea about the definition of these new subsets is that it is 
interesting to study not only the points where the resolvent of an operator is not de- 
fined, i.e. its spectrum, but also where this resolvent is large in norm. This explains 
the following definition of the e-pseudospectrum a e (A) of a matrix or an operator A, 



for any e > 0, if we write by convention that \\{zl — A) = +00 for every point z 
belonging to the spectrum a(A) of the operator. 

Let us mention that there exists an abundant literature about this notion of pseu- 
dospectrum. We refer here for the definition and some general properties of pseu- 
dospectra to the paper [15] of L.N. Trefethen. Let us also point out the more recently 
published book [16], which draws up a wide all-round view of this topic and gives a 
lot of illustrations. 

According to the previous definition, studying the pseudospectra of an operator is 
exactly studying the level lines of the norm of its resolvent. What is interesting in 
studying such level lines is that it gives some information about the spectral stability 
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of the operator. Indeed, pseudospectra can be defined in an equivalent way in term of 
spectra of perturbations of the operator. For instance, we have for any A € M n (C), 

a £ (A) ={z£C, ze a(A + B) for some B e M„(C) with ||S|| < e}. 

It follows that a complex number z belongs to the e-pseudospectrum of a matrix A if 
and only if it belongs to the spectrum of one of its perturbations A + B with ||_B|| < e. 
More generally, if A is a closed unbounded linear operator with a dense domain on a 
complex Hilbert space H, the result of Roch and Silbermann in [13] gives that 

a £ (A)= |J a(A + B), 

BeC(H), \\B\\ cm <s 

where C(H) stands for the set of bounded linear operators on H. From this second 
description, we understand the interest in studying such subsets if we want for example 
to compute numerically some eigenvalues of an operator. Indeed, we start to do 
it by discretizing this operator. This discretization and inevitable round-off errors 
will generate some perturbations of the initial operator. Eventually, algorithms for 
eigenvalues computing will determine the eigenvalues of a perturbation of the initial 
operator, i.e. a value in a e-pseudospectrum of the initial operator but not necessarily 
a spectral one. This explains why it is important in such numerical computations to 
understand if the e-pseudospectra of studied operators contain more or less deeply 
their spectra. 

Let us first notice that this study is a priori non-trivial only for non-self adjoint 
operators, or more precisely for non-normal operators. Indeed, we have for a normal 
operator A an exact expression of the norm of its resolvent given by the following 
classical formula (see for example (V.3.31) in [8]), 

(1.1.1) Vz^(A), \\{zl — A)" 1 1| = . 1 

d(z,a(A)) 

where d(z, c(A)) stands for the distance between z and the spectrum of the operator, 
when A is a closed unbounded linear operator with a dense domain on a complex 
Hilbert space. This formula proves that the resolvent of a normal operator cannot 
blow up far from its spectrum. It ensures the stability of its spectrum under small 
perturbations because the e-pseudospectrum is exactly equal in this case to the e- 
neighbourhood of the spectrum 

(1.1.2) <T £ {A) = {zEC:d(z,o-(A))<e}. 

Nevertheless it is well-known that this formula (1.1.1) is no more true for non-normal 
operators. For such operators, it can occur that their resolvents are very large in 
norm far from their spectra. This induces that the spectra of these operators can be 
very unstable under small perturbations. To illustrate this fact, let us consider the 
case of the rotated harmonic oscillator and the following numerical computation of its 
spectrum. The rotated harmonic oscillator is a simple example of elliptic quadratic 
differential operator 

H c = D 2 x + cx 2 , D x = r x d x , 
with c = e l7T / 4 . The numerical computation is performed on the matrix discretization 

((H c 'ifi, 5 , j)l 2 (R)) l<i,j<N' 

where N is an integer taken equal to 100 and (^j)j^m* stands for the basis of L 2 (R) 
composed by Hermite functions. The black dots appearing on this computation stand 
for the numerically computed eigenvalues. We can notice on this numerical simulation 
that the computed low energies are very close to theoretical ones since the spectrum 
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Figure 1 . Computation of some level lines of the norm of the resol- 
vent || (H c — = e^ 1 for the rotated harmonic oscillator H c with 
c = e 47r / 4 . The right column gives the corresponding values of log 10 s. 
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of the rotated harmonic oscillator is only composed of eigenvalues regularly spaced 
out on the half-line e iw / s R* + , 

a{H c ) = {e" /8 (2n + 1) : n G N}. 

However we notice that it is no more true for the high energies. It occurs for them 
some strong spectral instabilities, which lead to the computation of "false eigenvalues" 
far from the half- line e i7r / 8 lR?j_. Let us mention that some comparable computations 
can be found in [3]. In this paper, we are interested in studying when and how this 
kind of phenomena occurs in the class of elliptic quadratic differential operators. 

1.2. Elliptic quadratic differential operators. We study here the class of elliptic 
quadratic differential operators. It is the class of pseudodifferential operators defined 
in the Weyl quantization 

(1.2.1) q(x,O w u(x) = [ e ^-y^ q (^,Au(y)dyd^ 

by some symbols q(x, £), where (x, £) G M n x W 1 and n £ W, which are some complex- 
valued elliptic quadratic forms i.e. complex- valued quadratic forms verifying 

(1.2.2) (x, G M" x R", q(x, £) = (x, £) = (0, 0). 

Let us first notice that since the symbols of these operators are some quadratic forms, 
these arc only some differential operators, which are a priori non-self adjoint because 
their Weyl symbols are complex-valued. As mentioned before, the rotated harmonic 
oscillator is an example of such an operator since we have 

D 2 X + e w x 2 = (£ 2 + e t9 x 2 ) w , < 9 < tt, 

if D x = i~ 1 d x . This operator is a very simple example of non-selfadjoint operator for 
which we have noticed on the previous numerical simulation that it occurs some strong 



4 



spectral instabilities under small perturbations for its high energies. These phenomena 
have been studied in several recent works. We can mention in particular the works of 
L.S. Boulton [1], E.B. Davics [3], K. Pravda-Starov [10] and M. Zworski [18], which 
have given a good understanding of these phenomena. 

A question, which has been at the origin of this work, has been to study if these 
phenomena peculiar to the rotated harmonic oscillator are representative, or not, of 
what occurs more generally in the class of elliptic quadratic differential operators in 
every dimension. We have tried to answer to the following questions: 

- Does it always occur some strong spectral instabilities under small perturba- 
tions for the high energies of these operators ? 

- If it is not the case, is it possible to give a necessary and sufficient condition on 
the Weyl symbols of these operators, which ensures their spectral stability ? 

- Can we precisely describe the geometry, which separates the regions of the 
resolvent sets where the resolvents of these operators blow up in norm from 
the ones where one keeps a control on their sizes ? 

To understand these spectral stability or instability phenomena, we need to study 
the microlocal properties, which rule these phenomena in the class of elliptic quadratic 
differential operators. Let us mention that it is M. Zworski who first underlined in [18] 
the close link between these questions of spectral instabilities and some results of 
microlocal analysis about the solvability of pseudodifferential operators. 

1.3. Semiclassical pseudospectrum. To answer to these previous questions, it is 
interesting to use a semiclassical setting and to study a notion of pseudospectrum 
in this new setting. We define for a semiclassical family (Ph)o<h<i of operators on 
L 2 (R n ), with a domain D, the following notions of semiclassical pseudospectra. 

Definition 1.3.1. For all \i > 0, the set 

A s ;(A) = {z e C : VC > 0,V/i > 0,3 < h < h , \\{P h - z^W > CT^ 1 }, 

is called semiclassical pseudospectrum of index [i of the semiclassical family (Ph)o<h<i- 
The semiclassical pseudospectrum of infinite index is defined by 

With this definition, the points in the complement of the semiclassical pseudospectrum 
of index fi are the points of the complex plane where we have the following control of 
the resolvent's norm for sufficiently small values of the semiclassical parameter h, 

(1.3.1) 3C > 0,3h > 0,V < h < h , < eh-*. 

To prove the existence of semiclassical pseudospectrum of index /x, we will study the 
question of existence of semiclassical quasimodes 

(1.3.2) VC > 0, Vh a > 0, 3 < h< h , 3u h e D, 

\\uh\\mn") = 1 and \\P h u h - zu h \\ L 2 {Kn) < Ch^ 7 

in some points z of the resolvent set, which can be considered as some "almost eigen- 
values" in 0(h^) in the semiclassical limit. Let us notice that the definition chosen 
here for the notions of semiclassical pseudospectra differ from the one given in [5] for 
a semiclassical pseudodifferential operator. In fact, we have chosen a definition for 
semiclassical pseudospectra inspired by the remark made p. 388 in [5], because this 
definition only depends on the properties of the semiclassical operator rather than on 
its symbol. 
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The interest of working in a scmiclassical setting is a matter of geometry. We can 
explain this choice by the fact that it is easier for an elliptic quadratic differential oper- 
ator q(x, £) w to describe the geometry of semiclassical pseudospectra of its associated 
scmiclassical operator (q(x, h£) w )o<h<i, than to describe directly the geometry of its 
e-pseudospectra. The semiclassical setting is particularly well-adapted for the study 
of elliptic quadratic differential operators because there exists a simple link between 
this semiclassical setting and the quantum one. Indeed, using that the symbols of 
these operators are some quadratic forms q, wc obtain from the change of variables, 
y = h x ^ 2 x with h > 0, the following identity between the quantum operator q(x,£) w 
and its associated scmiclassical operator (q(x, h£) w )o < h<i, 

(1.3.3) q(x,Q w -l = ±(q(y,hr,) w -z), 

if z € C. This identity allows to get some information about the resolvent's norm 
behaviour of the quantum operator 

(q(x,O w -z)-\ 

if we have some information about semiclassical pseudospectra for its associated semi- 
classical operator. Let us mention for example that if a non-zero complex number z 
belongs to the scmiclassical pseudospectrum of infinite index of the operator 

(q(x,hO w )o<h<i, 

the identity (1.3.3) induces that the resolvent's norm of the quantum operator blows 
up along the half-line zK + with a rate faster than any polynomials 

(1.3.4) VJVeN,VC>0,V»jo>l,3»7>»j , \\(q(x,£) w - zrj)' 1 \\ > Cr] N , 

and this, even if this half-line zR + does not intersect the spectrum of the opera- 
tor q(x, £) w . Conversely, in the case where z £ A^ c (g(y, hrj) w ^j, z ^ and < \i < 1, 
the identity (1.3.3) shows that we can find some positive constants C\ and Ci such 
that the resolvent of the operator q(x, £,) w remains bounded in norm in some regions 
of the resolvent set of the shape 

(1.3.5) {ueC:\u\>C 1 ,d(A,u)<C 2 \proi A u\ 1 -"}n C\a(q(x,£) w ), 

where A = zM. + and proj A u stands for the orthogonal projection of u on the closed 
half-line A. Indeed, we obtain from (1.3.1) and (1.3.3) that 

3C> 0,3/70 >l,Vr?> 770, || (q(x,Q w - rje^)' 1 ]] < C rf~\ 

which induces that for all v e D(q{x, £) w ) and r/ > 770, 

\\{q(x,O w -Ve iaisz )v\\ LHRn) >C-y- M HU'(R»), 

if D(q(x,(;) w ) stands for the domain of the operator q(x,^) w . Then, we can find a 
constant 770 > 1 such that if z belongs to 

{ueC: \u\ > 770, d(e jarsz K + ,7i) < 2- 1 C , - 1 |proj eS „ B , R+ 7i| 1 -^} n C\a(q(x^) w ) , 
then 

|proj e ,„ g , K+ z| > 770. 

This induces using the previous estimates and the triangular inequality that if z 
belongs to 

{u e C : |u| > 770, d(e^ z R+,u) < 2- 1 C- 1 |proj el „ B , R+ 7 J | 1 -^} n C\a(q(x,O w ) , 
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we have for all v G D(q(x, S,) w ) , 

|| (q(x, W ~z)v\\ L2 > || (q(x, W - proj e ^ R+ ~z)v\\ L2 - d(e ia ^R + , z) \\v\\ L , 

> 2- 1 C- 1 |proj e » If! , R+ 5| 1 ^|M| i2 

> 2- 1 C- 1 t 1 1 ^\\v\\l^ 

because [i < 1. This last estimate shows that the resolvent of the operator q(x,£) w is 
bounded in norm by 2Ct]q~ 1 on the set 

{u G C : \u\ > fjo, d(e l ^ z R+,u) < 2- 1 C- 1 |proj e » Ig , R+U | 1 -^} n C\a(q{x,£) w ) . 

We notice that depending directly on the value of the index fi, < /j, < 1, the previous 
set contains more or less deeply in its interior the half-line 

{u G C : \u\ > fj , u G zR+}. 

This fact explains why in the following we will precise carefully the index of the 
semiclassical pseudospectrum to which a point does not belong when there is no 
semiclassical pseudospectrum of infinite index in that point. 

2. Statement of the results 

2.1. Some notations and some preliminary facts about elliptic quadratic 
differential operators. Let us begin by giving some notations and recalling known 
results about elliptic quadratic differential operators. Let q be a complex-valued 
elliptic quadratic form 

q : m; 1 x R£ -» C 

with n G N*, i.e. a complex- valued quadratic form verifying (1.2.2). The numerical 
range T,{q) of q is defined by the subset in the complex plane of all values taken by 
this symbol 

(2.1.1) S(g) = g(K£ x R£), 

and the Hamilton map F G M 2n (C) associated to the quadratic form q is uniquely 
defined by the identity 

(2.1.2) q((x, 0; (y, v)) = 0,F(y, V )) , (x, e R 2n , (y, V ) G K 2 ", 

where ?(•;•) stands for the polar form associated to the quadratic form q and a is the 
symplectic form on K 2 ™, 

(2.1.3) a((x, 0, (tf, »?)) = - ar.r?, 0, G K 2 ", (y, r?) G K 2 ". 

Let us first notice that this Hamilton map F is skew-symmetric with respect to a. 
This is just a consequence of the properties of skew-symmetry of the symplectic form 
and symmetry of the polar form 

(2.1.4) MX, Y G K 2n , a(X, FY) = q(X; Y) = q(Y; X) = a(Y, FX) = -a(FX, Y). 

Under this assumption of cllipticity, the numerical range of a quadratic form can 
only take some very particular shapes. It is a consequence of the following result 
proved by J. Sjostrand (Lemma 3.1 in [14]), 

Proposition 2.1.1. Let q : K™ x R|? — > C a complex-valued elliptic quadratic form. 
If n > 2, then there exists z G C* such that Re(zg) is a positive definite quadratic 
form. If n — 1, the same result is fulfilled if we assume besides that £(<?) ^ C. 
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This proposition shows that the numerical range of an elliptic quadratic form can only 
take two shapes. The hrst possible shape is when E(g) is equal to the whole complex 
plane. This case can only occur in dimension n = 1. The second possible shape is 
when £(<?) is equal to a closed angular sector with a top in and an opening strictly 
lower than ir. 

Figure 2. Shape of the numerical range T,(q) when S(g) ^ C. 
A 




Indeed, if T,(q) ^ C, using that the set T,(q) is a semi-cone 

tq(x, = q(Vtx, y/i$), t G K+, {x, f ) G R 2 ™, 
because g is a quadratic form, we have 

S(g) = M+z" 1 /, 

if z is the non-zero complex number given by the proposition 2.1.1 and / is the compact 
interval 

I=l + i lm(zq)(K), 
where K is the following compact subset of R 2n , 

|(x,e)eR 2 ":Re(zg)(x,a = l}. 

The compactness of if is a direct consequence of the fact that Re(zq) is a positive 
definite quadratic form. 

Elliptic quadratic differential operators define some Fredholm operators (see Lemma 3.1 
in [6] or Theorem 3.5 in [14]), 

(2.1.5) q(x,O w ■ L 2 (R"), 
where B is the Hilbert space 

(2.1.6) {u e L 2 (R n ) : x a D^u e L 2 (R n ) if \a + 0\ < 2}, 
with the norm 

\a+/3\<2 

The Fredholm index of the operator q(x, £) w + z is independent of z and is equal to 
if n > 2. In the case where n = 1, this index can take the values —2, or 2. More 
precisely, this index is always equal to if E(g) ^ C. 
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In the following, we will always assume that 7^ C. Under this assumption, 

J. Sjostrand has proved in the theorem 3.5 in [14] (see also Lemma 3.2 and Theorem 3.3 
in [6]) that the spectrum of an elliptic quadratic differential operator 

q(x,O w :B^L 2 (R n ), 

is only composed of eigenvalues with finite multiplicity 

(2.1.7) a(q(x,O w )={ J2 (rx + 2k x )(-i\):k x efi}, 

Aecr(F), 
- l AGS(g)\{0} 

where F is the Hamilton map associated to the quadratic form q and r\ is the dimen- 
sion of the space of generalized eigenvectors of F in C 2 ™ belonging to the eigenvalue 
A G C. Let us notice that the spectra of these operators is always included in the 
numerical range of their Weyl symbols. 

To end this review of preliminary properties of elliptic quadratic differential oper- 
ators, let us underline that the property of normality in this class of operators can be 
easily checked by computing the Poisson bracket of the real part and the imaginary 
part of their symbols 

<9Re q dim q dRc q dim q 

(2.1.8) {Req,lmq} = —.— ^— gp- 



Proposition 2.1.2. An elliptic quadratic differential operator 

q(x,t) w : B -> L 2 (R"), n G N*, 

is normal if and only if the quadratic form defined by the Poisson bracket of the real 
part and the imaginary part of its symbol is equal to zero 

(2.1.9) V(ar, G M 2n , {Re q, Im q}(x, £) = 0. 

Proof of Proposition 2.1.2. This proposition is a direct consequence of the composition 
formula in Weyl calculus (see Theorem 18.5.4 in [7]), which induces that the Weyl 
symbol of the commutator 

[q w , (q w T) = [Q w ,r] = -2i[(Re q) w , (Im q) w ], 

is equal to 

— 2z(Rc q $ Im q — Im q # Rc q) = — 2{Rc q, Im q}, 

because Re q and Im q are some quadratic forms. The notation Re q jj Im q stands 
for the Weyl symbol of the operator obtained by composition (Req) w (lmq) w . □ 

Remark. Let us notice that the symplectic invariance of the Poisson bracket (see 
(21.1.4) in [7]), 

(2.1.10) {(Re q) o x , (Im q) o x } = {Re q, Im q} o x , 

if x stands for a linear symplectic transformation of R 2n , implies that the condition 
(2.1.9) is symplectically invariant. 
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2.2. Statement of the main results. Let us consider an elliptic quadratic differ- 
ential operator 

q(x,O w :B^L 2 (R n ). 

We know from (2.1.7) that the spectrum of this operator is contained in the numerical 
range of its symbol The following proposition gives a first localization of the 

regions where the resolvent can blow up in norm and where spectral instabilities can 
occur. 

Proposition 2.2.1. Let q : R n x R n — > C, n e N*, be a complex-valued elliptic 
quadratic form. We have 

£<«,), || W ,, ._ l )-||<_ f _| j5J , 
where d(z, £(<?)) stands for the distance from z to the numerical range £(#). 

This result shows that the resolvent of an elliptic quadratic differential operator 
cannot blow up in norm far from the numerical range of its symbol. We are now 
going to study what kind of phenomena can occur in this particular set. There are 
two cases to separate according to the property of normality or non-normality of the 
operator. 

2.2.1. Case of a normal operator. Let us consider a normal elliptic quadratic differ- 
ential operator 

q{x,£) w : B -» L 2 (R n ). 

Let us recall that according to the proposition 2.1.2 this property of normality is 
exactly equivalent to the fact that 

V(z,0 e R 2n , {Re q,lm q}(x,£) = 0. 

In this case, we have the classical formula (1.1.1) for its resolvent's norm 

(2.2.1) V^M^D, -'>"' II- ().)) ■ 

which induces that the e-pseudospectrum of this operator is exactly equal to the 
e-neighbourhood of its spectrum 

W ) = C : d(z,a(q(x,O w )) < e}, e > 0. 

This classical formula (2.2.1) ensures that the resolvent cannot blow up in norm far 
from the spectrum and induces that the spectrum of such an operator is stable under 
small perturbations. 

Example 1. The operator 

(2.2.2) qi {x, W = -(1 + i)d 2 Xl - d 2 X2 + 4(-l + i) Xl 8 Xl + 2(-l + i)x 2 d Xl + 6ix 2 d X2 

+ 2ix l d X2 + (6 + hi)x\ + (11 + i)x\ + (10 + 4i)xix 2 -2 + 5i, 

is an example of a normal elliptic quadratic differential operator. Its spectrum is given 
by 

o-(qi(x,O w ) = {{2k! + 1) + (2k 2 + l)V2e^ : (k u k 2 ) G N 2 }. 
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FIGURE 3. Spectrum and a e-pseudospectrum of the operator qi(x,£) w . 
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Example. 2. Let us notice that when the numerical range is reduced to a closed 
half-line, the elliptic quadratic differential operator q(x,£) w is always normal since 

{Re q, Im q} = |z| 2 {Rc(z" 1 g), ^(z" 1 ^)} = 0, 

if z G C* is chosen such that Im(z~ 1 q) = 0. In fact, the operator q(x,^) w can in this 
particular case be reduced after a conjugation by a unitary operator on L 2 (R") to the 
operator 

n 
3=1 

where Xj > for all j = 1, n. 

FIGURE 4. Example of a normal elliptic quadratic differential operator. 
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2.2.2. Case of a non-normal operator. Let us consider a non-normal elliptic quadratic 
differential operator 

q(x,O w :B^L 2 (R n ), nef. 

We assume in the following that the numerical range £((/) is distinct from the whole 
complex plane 

(2.2.3) £(<?)^C. 

As mentioned in the section 2.1, this additional assumption is always fulfilled in 
dimension n > 2. It only excludes some very particular one-dimensional elliptic 
quadratic differential operators (see the remark following the proposition 2.2.2 for 
more precision about these operators). 

Under this additional assumption, the numerical range £(<?) is always a closed 
angular sector with a top in and a positive opening strictly lower than tt. 

2.2. 2. a. On the pseudospectrum at the interior of the numerical range. Let us consider 
the associated semiclassical elliptic quadratic differential operator 

{q(x,hO w )o<h<i. 

We can build in every point of the interior of the numerical range Ti(q) some semi- 
classical quasimodes. 

Theorem 2.2.1. If the elliptic quadratic differential operator 

q(x,£) w :B^L 2 {R n ), n E N* , 

is non- normal and verifies £(<?) ^ C then for all z G £(g) and JVel, there exist 
ho > and a semiclassical family (uh)o<h<h € S(M. n ) such that 

ll^hlU^R") = 1 and \\q(x, h^) w u h - zu h \\ L 2 {Rn) = 0(h N ) when h -> + . 

This result induces the existence of semiclassical pseudospectrum of infinite index in 
every point of the interior of the numerical range E(g). 

According to (1.3.4), this result in the semiclassical setting induces that the resol- 
vent's norm of the quantum operator q(x,£,) w blows up fastly along all the half-lines 
belonging to the interior of the numerical range £(<?), 

(2.2.4) V^eE(g),VJVeN,VC>0,V%>l,3»7>»jd, \\(q(x,Z) w - zrf)' 1 ]] > Cr] N . 

We deduce from (2.1.7) that as soon as an elliptic quadratic differential operator is 
non-normal its resolvent blows up in norm in some regions of the resolvent set far 
from its spectrum. This fact induces that the high energies of such an operator are 
very unstable under small perturbations as we have already noticed on the numerical 
computation performed for the rotated harmonic oscillator. It follows that in the class 
of elliptic quadratic differential operators 1 the property of spectral stability is exactly 
equivalent to the property of normality: 

cr(q(x,£) w ) is stable under q (x, £) w is a normal 4^ {Rc q, Im g} = 0. 

small perturbations operator 

By spectral stability, we mean here that the resolvent of these operators cannot blow 
up in norm far from their spectra. Let us add that it is not very surprising to have 
this property of spectral stability under the assumption of normality, but it is worth 



If we exclude the one-dimensional particular cases previously mentioned. 
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noticing that as soon as this property is violated, it occurs in this class of operators 
some strong spectral instabilities under small perturbations for their high energies. 

Examples. The two following operators 

(2.2.5) q 2 (x,Z) w = -d 2 Xl -2d 2 X2 +Aix 2 d X2 +2x\ + {A + i)x 2 2 + Ax l x 2 + 2i 
and 

(2.2.6) q 3 (x, W = -(1 + i)d 2 Xl - 2d 2 X2 + 4(-l + i) Xl d Xl + 2(1 - i)x 2 d Xl - Ai Xl 8 X2 

+ (9 + Ai)x\ + (2 + i)x\ - 4(1 + i)x x x 2 -2 + 2i, 
are some examples of non-normal elliptic quadratic differential operators. 

2.2.2.b. On the pseudospectrum at the boundary of the numerical range. Let us now 
study what occurs on the boundary of the numerical range dT,(q) for a non-normal 
elliptic quadratic differential operator 

q(x,£) w : B -» L 2 (R"). 

Let us mention that we always assume that ^ C. Under these assumptions, the 
boundary of the numerical range is composed of the union of the origin and two 
half-lines Ai and A 2 , 

(2.2.7) 0£(g) = {0}U A 1 U A 2 , 
that we can write 

(2.2.8) Ai = ZiR* + and A 2 = z 2 R* + with z u z 2 G dT,(q) \ {0}. 

We need to define a notion of order for the symbol q(x, £) on these two half- lines Aj, 
j = 1, 2. Let us begin by recalling the classical definition of the order fc(xo,^o) of a 
symbol p(x, £) at a point (a; ,Co) € IR 2 ™ (see section 27.2, chapter 27 in [7]). This 
order k(x , £ ) is an element of the set N U {+00} defined by 

(2.2.9) k(x ,t ) = sup{3 eZ:pi(x ,t ) = 0, Vl<|7|<j}, 

where I — i 2 , ik) <E {l,2} fe , |7| = k and pi stands for the iterated Poisson 
brackets 

Pi = H Pii Hp i2 ...Hp iki pi k , 

where p\ and p 2 are respectively the real and the imaginary part of the symbol p, 
p = Pi + ip 2 . The order of a symbol q at a point z is then defined as the maximal 
order of the symbol p = q — z at every point (xo, £0) S M 2 ™ verifying 

p(xo,£,o) = q(xo,£,o) -z = Q. 

Let us underline that the symplectic invariance of the Poisson bracket (2.1.10) induces 
the same property for the order of a symbol at a point. 

Since here the symbol q is a quadratic form, all the iterated Poisson brackets are 
also some quadratic forms. This property of degree two homogeneity of these Poisson 
brackets induces that the symbol q has the same order at every point of each half-line 
Aj, j = 1,2. This allows to define the order of the symbol q on the half- line Aj by 
defining this order by this common value. Let us mention that this order can be finite 
or infinite. 

Examples. One can easily check that the Weyl symbol 

i 2 + e l0 x 2 , O<0<tt, 
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of the rotated harmonic oscillator has an order equal to 2 on the both half-lines R!j_ 
and e l6 R+, which composes the boundary of its numerical range. The symbol (72 of 
the operator defined in (2.2.5) has an order equal to 2 on iR!j_ and to 6 on , 

E(qr 2 ) = {z e C : Re z > 0, Im z > 0}. 

On the other hand, we can verify that the symbol q$ of the operator defined in (2.2.6) 
is of infinite order on the half-line R+ and has an order equal to 2 on e 47r//4 M!j_, 

s (3s) = {0} U {z e C* : < arg z < tt/4}. 

In the case where the symbol is of finite order on a half-line Aj, j = 1, 2, we have 
the following result. 

Theorem 2.2.2. If the Weyl symbol q(x , £) of a non-normal elliptic quadratic differ- 
ential operator is of finite order kj on the half-line 

Aj, je{l,2}, Aj c dT,(q) \ {0}, 

then this order is necessary even and there is no semiclassical pseudospectrum of 
index kj/(kj + 1) on Aj for the associated semiclassical operator 

A J cC\A- /(fej+1) (^,^r)- 

Remark. Let us mention that we can more precisely establish that in dimension n > 1, 
the order kj is an even integer verifying 

2 < kj < An - 2. 

This result is proved in [12]. 

By rephrasing this result in a quantum setting, it follows from (1.3.5) and (2.1.7) 
that when the symbol q of a non-normal elliptic quadratic differential operator q(x, !;) w 
is of finite order kj on a half-line 

Aj, je{l,2}, Aj c dT,(q) \ {0}, 

then the resolvent of this operator remains bounded in norm in a set of the following 
type 

(2.2.10) {ue C : |it| > Ci, d{Aj,u) < C 2 |proj A .«| W}, 

where C\ and Ci are some positive constants. 

As we will see in its proof, this absence of semiclassical pseudospectrum is linked 
to some properties of subellipticity. Let us just underline for the moment that the 
index kj/(kj + 1), which appears in this result is exactly equal to the loss appearing 
in the subelliptic estimate hidden behind this result. 

About the case of infinite order, the situation is much more complicated. Never- 
theless, we can first notice in this case that we cannot expect to prove a stronger result 
than an absence of semiclassical pseudospectrum of index 1. Indeed, we can easily 
check on the example of the operator (73(21, £) w defined in (2.2.6) that its spectrum is 
given by 

a(q 3 (x,O w ) - {(2k! + 1)V2+ (2fc 2 + l)3^e*t : (k u k 2 ) € N 2 }. 

We recall that the spectrum of this operator is only composed of eigenvalues and that 
its symbol is of infinite order on . It follows from the structure of the spectrum and 
(1.3.5) that if there is no semiclassical pseudospectrum of infinite index in a point of 
the half-line R^, there is necessary no semiclassical pseudospectrum of index /j, with 
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an index fi > 1. In fact, we can prove by using a result of exponential decay in time 
for the norm of contraction semigroups generated by elliptic quadratic differential 
operators (see [12]) that there is never some semiclassical pseudospectrum of index 1 
on all these half-lines of infinite order. Let us mention that this result of exponential 
decay will not be proved here but it will be explained in the following how it induces 
the absence of semiclassical pseudospectrum of index 1. 

2.2.3. About the geometry of s-pseudospectra for elliptic quadratic differential opera- 
tors. Let us now explain what are the consequences of these results on the geometry 
of e-pseudospcctra for elliptic quadratic differential operators. Let us begin by con- 
sidering the one-dimensional case which is a bit particular. In dimension n = 1, an 
elliptic quadratic differential operator can be reduced after a similitude and a conju- 
gation by a unitary operator to the harmonic oscillator or to the rotated harmonic 
oscillator. 

Proposition 2.2.2. Let us consider q :Rxl^C a complex-valued elliptic quadratic 
form such that S(q) =/= C. For all h > 0, there exist a unitary operator (more precisely 
a metaplectic operator) Uh on L 2 (R), which is an automorphism of the spaces <S(R) 
and B, z £ C* and 9 £ [0, 7r[ such that 

Vh>0, q{xMT =zU h ((hD x f + e i6 x 2 )U^. 

Remark. In the case where = C, an elliptic quadratic differential operator 

q(x, £)™ can be reduced after a similitude and a conjugation by a unitary operator on 
L 2 (R n ) to the operator defined in the Weyl quantization by the symbol 

(£ + ix)(£ + nx) with n £ C, Im 77 > 0, 

or 

(£ — ix)(t; + nx) with n £ C, Im n < 0, 
depending on the value of its Fredholm index, which is equal to —2 in the first case 
and to 2 in the second one. 

As we will see in the following, this proposition allows us to reduce the study of a 
one-dimensional non-normal elliptic quadratic differential operator verifying 

£(«) + C, 

to the one of the rotated harmonic oscillator 

H e = D 2 x + e l6 x 2 , < 6 < ir. 

Let us mention that the previous results (Theorem 2.2.1 and Theorem 2.2.2) were 
already known in the particular case of the rotated harmonic oscillator. Indeed, 
the existence of semiclassical quasimodes inducing the presence of semiclassical pseu- 
dospectrum of infinite index in every point of the interior of the numerical range for 
the associated semiclassical operator, is a direct consequence of a result proved by 
E.B. Davies in [4] (Theorem 1). About the absence of semiclassical pseudospectrum 
of index 2/3 on the boundary of the numerical range, this result has been proved for 
the rotated harmonic oscillator in [10] 2 . 

As proved in [10], this absence of semiclassical pseudospectrum allows to give a 
proof of a conjecture stated by L.S. Boulton in [1]. It deals with the geometry of e- 
pseudospcctra for the rotated harmonic oscillator. Let us now recall some facts about 
this conjecture and some results proved by L.S. Boulton in [1]. 



Let us recall that the value of the order is equal to 2 in this case. 
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L.S. Boulton has first proved (Theorem 3.3 in [1]) that the resolvent of the rotated 
harmonic oscillator blows up in norm along all a family of curves of the following form 

■q ^ b-q + e ie if , 

where b and p are some positive constants verifying 1/3 < p < 3, 

(2.2.11) \\(H e - {br) + e ie rj p ))~ 1 \\ -> +oo when r, -» +oo. 

On the other hand, he also proved that the resolvent of this operator remains bounded 
in norm on two half-stripes parallel to the half- lines R + or e l6 M. + . More precisely, he 
proved that there exist some positive constants d and M d such that 

(2.2.12) sup \\{H -{7 1 + ib)Y 1 \\<M d , 

(2.2.13) sup -e w (j?-*&)) _1 || < M d . 

These bounds provide some information about the shape of e-pseudospectra of the 
operator Hg. Indeed, L.S. Boulton has proved using these results that for all suffi- 
ciently small value of the positive parameter e, the e-pseudospectra of the rotated 
harmonic oscillator is contained in the shaded set appearing on the following figure. 
The eigenvalues appear on this figure marked by some o. 

Figure 5. A first localization of the e-pseudospectra of the rotated 
harmonic oscillator. 




More precisely, L.S. Boulton proved that for all < S < 1 and m E N, there exists 
a positive constant eo such that for all < e < eo, 

m 

(2.2.14) a e (H g ) c\J{ze C:\z- X n \ <S}U [A m+1 - 6e w ' 2 + Sg] , 

n=0 

where 

A„ = e 4e / 2 (2n + l), neN 

and 

Sg = {z e C* : < arg z < 9} U {0}. 
In fact, in view of some numerical calculations performed by E.B. Davies in [3], 
L.S. Boulton has conjectured that the index p = 1/3 appearing in (2.2.11) is the 
critical one in the following sense: 
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Let us consider < p < 1/3, < 5 < 1 and m e N. If b m , p and E are some positive 
constants verifying 

b m .,pE + e l6 E p = X m and Vr/ > E, arg ^ < (9/2, 

where z,, = b m ^n + e l0 -q p , let us set 

tt m ,p = {\z v \e la € C : »y > £, arg ^ < a < arg(z^e l9 )}. 

L.S. Boulton has conjectured the following result. 

Boulton's conjecture. There exists e > such that for all < e < e 0l 

m 

(2.2.15) a £ {Hg) C (J {z G C : |z - A„| < 5} U n m , p . 

n=0 

The absence of semiclassical pseudospectrum of index 2/3 on the boundary of the 
numerical range dT,(q) \ {0} for the rotated harmonic oscillator 3 given by the theorem 
2.2.2 shows that this index 1/3 is actually the critical one. Indeed, we can deduce 
(2.2.15) from (2.2.10) (sec [10] for more details) since here % = 2, j e {1,2}. As 
we will see, this theorem 2.2.2 is a consequence of a subelliptic estimate for gen- 
eral semiclassical pseudodifferential operators proved by N. Dencker, J. Sjostrand and 
M. Zworski in [5] (Theorem 1.4). In the particular case of the rotated harmonic oscil- 
lator, a more elementary proof of this result using only some non-trivial localization 
scheme in the frequency variable is given in [10]. 

Let us notice that this inclusion (2.2.15) allows to give a sharp description of the e- 
pseudospectra of the rotated harmonic oscillator, which is optimal in view of (2.2.11). 



Figure 6. Shape of the e-pseudospectra of the rotated harmonic oscillator. 




By coming back to the case of an arbitrary dimension n > 1, let us finally underline 
that using the theorem 2.2.2, we can give similar descriptions of the £-pseudospectra 
for non-normal elliptic quadratic differential operators, to the one given by L.S. Boul- 
ton for the rotated harmonic oscillator, when the symbols of these operators are of 
finite order on the two open half-lines, which compose the boundary of their numerical 
ranges. The only difference with the particular case of the rotated harmonic oscillator 
is that the critical indices, which appear in this description can be different. Indeed, 



'The order of the rotated harmonic oscillator's symbol is equal to 2 on dH(q) \ {0}. 
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these critical indices depend directly according to (2.2.10) on the order of the symbols 
on the two half-lines composing the boundary of their numerical ranges. We refer the 
reader to [10] for more details about the way of getting from (2.2.10) such descriptions 
of e-pseudospectra. 

3. The proofs of the results 

Before giving the proofs of the results stated in the previous section, let us begin by 
recalling the symplectic invariance property of the Weyl quantization (see Theorem 
18.5.9 in [7]). This symplectic invariance is actually the most important property of 
the Weyl quantization. 

For every affine symplectic transformation \ of R 2 ™, there exists a unitary trans- 
formation U on L 2 (R n ), uniquely determined apart from a constant factor of modulus 
1, such that U is an automorphism of the spaces <S(R"), B and 6>'(R™), where B is 
the Hilbert space defined in (2.1.6), and 

(3.0.1) (ao X )(x,O w = U- 1 a(x,O w U, 

for all a e <S'(R 2 ™). The operator U is a metaplectic operator associated to the affine 
symplectic transformation \- 

This symplectic invariance of the Weyl quantization induces the same property for 
the semiclassical pseudospectra of elliptic quadratic differential operators in the sense 
that if 

q : R™ x R™ -» C, 

is a complex- valued elliptic quadratic form and x is a linear symplectic transformation 
of R 2 ™, we have for all /i G [0, oo], 

(3.0.2) A^((qo X )(x,hO w )=A^(q(x,H) w ). 

To prove this fact, let us begin by noticing that for all a e <S'(R 2n ) and ft, > 0, we 
have 

U^a(x,O w U h = aih-^x^tT, 

where 

U h f{x) = h^fih^x), 

since according to the proof of Theorem 18.5.9 in [7], Uh is a metaplectic operator 
associated to the linear symplectic transformation 

(x,0^(h- 1/2 x,h^ 2 0- 

Let us now consider the case where the symbol a is a quadratic form. The homogeneity 
property of such a symbol implies that 

V/i > 0, a(h-^ 2 x, h 1 ' 2 ^) = ^-a(x, h£), 

and 

Vft > 0, U^a{x,O w U h = ^a{xMT- 

If q : R™ x R^ — > C is a complex-valued elliptic quadratic form and x is a linear 
symplectic transformation of R 2 ™, we can notice that 

(qo X )(x,hZ) w , h>0, 

is actually an elliptic quadratic differential operator since the symbol qo X is an elliptic 
quadratic form. Let zeC and U be a metaplectic operator associated to the linear 
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symplectic transformation \- Using that U and Uh are some automorphisms of the 
Hilbert space B and 

(3.0.3) U^U^UhqixMTU^UUh = U^U^hqix^yUUn 

= hU^(q o X )(x, $ w U h = (q o X )(x, hZ) w , 

we obtain that 

U^U-'U^qix, h£) w - zy'u^UUn = ((q o X )(x, h£) w - z)~\ 

Using finally that U^ 1 U^ 1 Uh is a unitary transformation of L 2 (R"), this identity 
implies that 

WHqox^htr-z)- 1 ]] = \\(q(x,hO W -zT'l 
which proves (3.0.2). In the following, this property of symplectic invariance will 
allow us to reduce certain symbols to some normal forms by choosing new symplectic 
coordinates. We can now begin to prove the results stated in the previous section. 

Let us start by the proof of the proposition 2.2.1. 

Proof of Proposition 2.2.1. If the numerical range is equal to the whole complex plane, 
there is nothing to prove. If S(q) ^ C, we have seen in the previous section that the 
numerical range is necessary a closed angular sector with a top in and an opening 
strictly lower than it. 

Let us consider z and denote by Zq its orthogonal projection on the non- 

empty closed convex set According to the shape of the numerical range, it 

follows that z belongs to its boundary and that we can find a complex number 
Z\ £ C*, \zi\ = 1 such that 

T,(z!q) C {z £ C : Re z > 0} 

and 

(3.0.4) ziz £ {z £ C : Re z < 0}, d(z, £(<?)) = d(z lZ , iR). 

Using now that the operator i[Im(zig)] 1 " is formally skew-selfadjoint, we obtain that 
for all u £ S(R n ), 

Re(ziq(x, £) w u - z-lzu, u) L2(Rn) 
(3.0.5) = d(^,«)|M|| 2(Rfl) + ([Re(z iq (x,Z))] w u,u) L2(Rny 

Then, since the quadratic form Re(ziq) is non-negative, we deduce from the symplectic 
invariance of the Weyl quantization and the theorem 21.5.3 in [7] that there exists a 
metaplectic operator U such that 

k k+l 

[Rc(z iq ^0)] W = r 1 (^A i (^+^)+ £ x^)u, 

j=i j=k+i 

with k, I £ N and Xj > for all j = 1, k. By using that U is a unitary operator on 
L 2 (R"), we obtain that the quantity 

([Rc(z iq {x,0)] W u,u) L2(Rn) 

k k+l 

= ^2 X A\\ D Xl Uu \\h(^) + \\ x 3 Uu \\h(M^)) + \\ X 3 Uu \\ 2 LHMn), 
3=1 j=k+l 
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is non-negative. Then, we can deduce from the Cauchy-Schwarz inequality, (3.0.4) 
and (3.0.5) that for all u e S(R n ), 

d{z,E(q))\\u\\ L 2 {R n ) < N \\q{x,£) w u- Zu\\ L 2( R ny 

Finally, using the density of the Schwartz space 5(R") in B and the fact that \z\\ = 1, 
we obtain that 

v,*e(,), n W ,,a.-,)-||< 5{ - ; i j5ri 

since according to (2.1.7), a(q(x, £) w ) C T,(q). □ 

We now consider the one-dimensional case, which is a bit particular. 

3.1. The one-dimensional case. In dimension n = 1, we can reduce the study of 
complex- valued elliptic quadratic forms to exactly three normal forms after a simili- 
tude and a real linear symplectic transformation. 

Lemma 3.1.1. Let q : R x x — > C be a complex-valued elliptic quadratic form in 
dimension 1 . Then, there exists a linear symplectic transformation \ of M 2 such that 
the symbol q o \ is equal to one of the following normal forms: 

(i) a(£ 2 + e i9 x 2 ) withaeC*, < (9 < tt. 

(ii) a(£ + ix)(£ + r]x) with a e C*, rj e C, Im 77 > 0. 
(in) a(£ — ix)(£ + tjx) with a e C* , i] e C, Im r\ < 0. 

In the two last cases (ii) and (Hi), the numerical range is equal to the whole 

complex plane, £(<?) = C. 

Proof of Lemma 3.1.1. Let q : M 2 — > C be a complex- valued elliptic quadratic form. 
Let us first consider the case where £(<?) 7^ C. We deduce from the proposition 2.1.1 
that we can reduce our study to the case where Re q is a positive definite quadratic 
form. Then, using Lemma 18.6.4 in [7], we can find a real linear symplectic transfor- 
mation to reduce the quadratic form Re q to the normal form 

X(x 2 +£ 2 ), with A > 0. 

It follows that there exist some real constants a, b and c such that 

q{x,£) = \(x 2 +C 2 +i(ax 2 + 2bx£ + c£, 2 )). 

Then, we can choose an orthogonal matrix P € 0(2, M.) diagonalizing the real sym- 
metric matrix associated to the quadratic form ax 2 + 2bx£ + c£ 2 , 

a b\ ( \i 
be) \ A 2 

with Ai,A 2 e E. If P € 0(2, R) \ 50(2, E), we have 

a b\~ = ( \ 2 
be) \ Ai 

if (To is the matrix with determinant equal to —1, 

1 

1 

and P = Pao. It follows that we can always diagonalizc the real symmetric matrix 
associated to the quadratic form A _1 Im q by conjugating it by an element of 50(2, R). 



20 



Since the symplectic group is equal in dimension 1 to the group SL(2, R), we can after 
a linear symplectic transformation of R 2 reduce the quadratic form q to 

X(x 2 +C 2 + i{ lx x 2 + 72 £ 2 )) = a(f + re ie x 2 ), 

where 71,72 € R, a € C*, r > and e] - 7r,7r[. Let us notice that the elliptic- 
ity of q actually implies that 9 ^ 7r[27r]. Finally, using the real linear symplectic 
transformation (a;,£) (r _1 / 4 x, r 1 / 4 ^), we get a symbol of type (i), 

if < ^ < 7r. If — tv < 8 < 0, we need to use besides the real linear symplectic 
transformation (a;, £) (£, — a;) to obtain a symbol of type (i), 

arh ie (e + e- ie x 2 ). 

Let us now assume that = C. Since the dimension is equal to 1, we can factor the 
symbol q on C as a polynomial function of degree 2 in the variable £. Thus, according 
to the dependence in the variable x of the polynomial function's coefficients, we can 
find some complex numbers Ai, A2 and a G C* such that 

q(x, £) = a(£ - A^)^ - A 2 x). 

The ellipticity assumption for the quadratic form q induces that 

Im Xj ^ 0, 

if j = 1, 2. Using now the linear symplectic transformation (x, £) 1— > (a;, £ + Re Aia;), 
we can assume that 

(3.1.1) gfoO^tf-iraOte + kr), 

with r € R* and Im 6 ^ 0. Let us now check that the assumption S(q) = C induces 
that r Im b < 0. Since 

(£ — irx){^ + 6a;) = £ 2 + (b — ir)x^ — irbx 2 , 

the condition S(g) = C implies that for all (v,w) e R 2 , there exists a solution 
(xo, Co) € R 2 of the system 

£ 2 + Re 6 x£ + r Im 6 a; 2 = v 
a;£(Im b — r) — r Re b x 2 = w. 

Let us first notice that the second equation of (3.1.2) is fulfilled for all w e R only if 

Im b r. 

If w ^ 0, it follows from the second equation of (3.1.2) that a; 7^ and 

(3.1.3) 4o = -7f — r x • 

(Im b — r)x 

Let us consider the case where v = 0. Using (3.1.3) and the first equation of (3.1.2), 
we obtain that 

(w + r Re b a; 2 ,) 2 + Re b (Im b - r)xl(w + r Re b x 2 Q ) + r Im b (Im b - rfx^ = 0. 
We can rewrite this equation as f w (Xo) = if we set Xq = Xq and 

(3.1.4) f w (X) =rlmb ((Re b) 2 + (Im b - r) 2 )X 2 + w Re b (Im b + r)X + w 2 . 

Thus, the condition £(g) = C implies that there exists for all w/0, a non-negative 
solution X of the equation f w (X ) = 0. Since the quantity r Im b is assumed to be 
non-zero, we first study the case where r Im b > 0. In this case, since 

(3.1.5) f w {X) = 2rlmb ((Re bf + (Im b - rf)X + w Re b (Im 6 + r) 



(3.1.2) 
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and 

2r Im b ((Re bf + (Im b - rf) > 0, 
because Im b ^ r, we have 

(3.1.6) VA e R+, / W (A) > / w (0) = W 2 > 0, 

if w 7^ and 

k; Re b (Im 6 + r) 

< 

2r Im 6 ((Re b) 2 + (Im b - r) 2 ) ~ 

The estimate (3.1.6) shows that if r Im b > 0, the equation f w (X) = has no non- 
negative solution for all value of the parameter w ^ 0. This proves that the condition 
£(<?) = C induces that r Im b < 0. Using the linear symplectic transformation 

{x,0»{\r\- 1/2 x,\r\ 1/2 0, 
we obtain the normal forms (ii) and (Hi), 

a\r\(t; + ix)(£ + i]x) with Im r\ > and a|r|(£ — ix)(£ + r]x) with Im 77 < 0, 

where 77 = |r| _1 &. Finally, we can easily check that the numerical ranges of the normal 
forms (ii) and (Hi) are actually equal to the whole complex plane C. □ 

Let us notice that the proposition 2.2.2 and the remark following its statement are 
some direct consequences of the symplectic invariance property of the Weyl quanti- 
zation (see (3.0.3)) and the previous lemma. We can add that as proved after the 
lemma 3.1 in [6], the Frcdholm indices of the one-dimensional elliptic quadratic dif- 
ferential operators with symbols of type (i), (ii) and (Hi) are respectively equal to 0, 
-2 and 2. 

As we have mentioned in the previous section, the results of Theorem 2.2.1 and 
Theorem 2.2.2 are already known in the particular case of the rotated harmonic oscil- 
lator. The existence of semiclassical quasimodes inducing the presence of semiclassical 
pseudospectrum of infinite index in every point of the interior of the numerical range 
for the associated semiclassical operator, is a direct consequence of a result proved 
by E.B. Davies in [4] (Theorem 1) and; the absence of semiclassical pseudospectrum 
of index 2/3 on the boundary of the numerical range has been proved for the ro- 
tated harmonic oscillator in [10] 4 . As we have previously mentioned (see (2.1.10) and 
(3.0.2)), the property of non-normality, the order of symbols and the semiclassical 
pseudospectra of elliptic quadratic differential operators are symplectically invariant. 
These properties allow us to reduce by any real linear symplectic transformations the 
symbols of the elliptic quadratic differential operators that we consider in our proof of 
the theorem 2.2.1 and the theorem 2.2.2. By using the lemma 3.1.1, we deduce from 
the results of the theorem 2.2.1 and the theorem 2.2.2 proved for the rotated harmonic 
oscillator that they are therefore also fulfilled by all non-normal one-dimensional el- 
liptic quadratic differential operators with a numerical range different from the whole 
complex plane. 

We now consider the multidimensional case. As we will see in the following, there is 
a real jump of complexity between the one-dimensional case and the multidimensional 
one. This jump is among other things a consequence of the complexity increase of 
symplectic geometry in dimension n > 2 and the larger diversity appearing in the 
class of elliptic quadratic differential operators. 



Let us recall that the value of the order is equal to 2 in this case. 
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3.2. Case of dimension n > 2. We only need to study the case of a non-normal 
elliptic quadratic differential operator 

(3.2.1) q(x,O w :B^L 2 (R n ), 

in dimension n > 2. Let us recall that in this case, the numerical range E(g) is a 
closed angular sector with a top in and a positive opening strictly lower than tt, and 
that the proposition 2.f .2 gives that 

(3.2.2) 3(3:0, &) G R 2 ", {Ro 9 , Im <z}(x , &) ^ 0. 

Let us begin by studying what occurs at the interior of the numerical range 

3.2.1. On i/ie pseudospectrum at the interior of the numerical range. To prove the 
existence of semiclassical quasimodes for the associated semiclassical operator given 
by the theorem 2.2.1, we need a first purely algebraic step to characterize the points 
belonging to the interior of the numerical range. 

Let us consider the following decomposition of the numerical range 

(3.2.3) E(q)=A\JB, 
where 

(3.2.4) A = {z e £(«) : 3(x ,£o) € R 2 ", z = q(x ,&), {Rc g,Im q}(x a ^ ) ? 0} 
and 

(3.2.5) B = {ze £(<?) : z = q(x ,^) => {Rc q,Im q}(x ,£ ) = 0}. 

The next section is devoted to give a geometrical description of these two sets. We 
establish using purely algebraic arguments that 

(3.2.6) I = ±(q) and B = dX(q). 

This result is a consequence of the geometry induced by the quadratic setting to which 
the studied symbols belong. 

Let us begin by noticing that the symplectic invariance of the Poisson bracket 
(2.1.10) induces the same property for the sets A and B. We can therefore use some 
real linear symplectic transformation to reduce the symbol q. Since 

{Re(zq), lm(zq)} = |z| 2 {Re q, Im q}, 

we deduce from this symplectic invariance, from the proposition 2.1.1 and the lemma 
18.6.4 in [7] that after a similitude, we can reduce our study to the case where 

n 

(3-2.7) Re q(x,0 =^A J -(^+^), 

with Xj > for all j = 1, ...,n. 

3. 2.1. a. Geometrical description of the sets A and B. We begin by proving the fol- 
lowing inclusion 

(3.2.8) 0E(g) C B. 

Let us consider z € dT,(q) and (x ,^o) G M 2 " such that z — q(xo,£o)- This is 
possible because the numerical range is a closed angular sector. If z = 0, the ellipticity 
property of q implies that 

(z ,£o) = (0,0) and {Re q,Im q}(x Q ,£, ) = 0, 

because this Poisson bracket is also a quadratic form. This proves that z € B. If 

z e flE(g) \ {0}, 
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let us consider the global solution Y of the linear Cauchy problem 

(3.2.9) jrw = /W^)) 



Oo,£o), 

associated to the Hamilton vector field of the symbol Re q, 



/<9Rc q d dRc q d 



d£j dxj dxj d£j 



)■ 



It is actually a linear Cauchy problem since Re q is a quadratic form. Setting 

f(t)=Imq(Y(t)), 

a direct computation gives that 

/'(0) = {Re«,Im q}(x ,£ ). 

If /'(0) ^ 0, we could find t ^ such that 

l/(*o)| > |/(0)| = |Im z\. 

Since Y is the flow associated to the Hamilton vector field of Re q, the quadratic form 
Rc q is constant under it. It follows that for all teR, 

Re q(Y(t)) = Re q(Y(0)) = Re z 

and provides a contradiction because, since z € d'E(q) \ {0}, this would imply in view 
of the shape of the numerical range S(g) (see Figure 7) that 

q(Y(t )) ? £(«). 

It follows that the Poisson bracket {Re q, Im <?}(xo,£o) is necessary equal to and 





Figure 7 


i 


i 




Q(Y(t )) . 
















that zeB. This ends the proof of the inclusion (3.2.8). 

Let us now assume that 

(3.2.10) 9£(g) C B, 0£(g) ^ B. 
In this case, we could find 

(3.2.11) zeB\dT,(q). 
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Let us first notice that z is necessary non-zero since € dT,(q), and that Re z > 0, 
since from (3.2.7), 

(3.2.12) %)\{0}c{zeC*:Rcz>0}. 
The fact that z belongs to the set B implies that 

(3.2.13) {ZtiZZl ^{^ q ,lm q}M = 0. 

We also know that there exists at least one solution to the system appearing in the 
left-hand-side of (3.2.13). Since from (3.2.7), the quadratic form Re q is positive 
definite, we can simultaneously reduce the quadratic forms Re q and Im q by finding 
an isomorphism P of K 2n such that in the new coordinates y — P^ 1 (x, £), 

2n In 

(3.2.14) Re q(Py) = ^ Vj and Im Q( p y) = ^2 a jVj with "i < - < 

Let us now consider the following quadratic form 

(3.2.15) p(y) = {Re g,Im q}(Py). 
We get from (3.2.13) and (3.2.14) that 

(3.2.16) |S l2/ ^ Re T Z =^p(y) = 0. 

Let us underline that the isomorphism P is not a priori a symplectic transformation 
and that it does not preserve the Poisson bracket {Re q,lm q}. 
We consider the two following sets 

(3.2.17) E 1 = {y e K 2 " : r(y) = 0}, 
where 

2n 

(3.2.18) r(y) = £ (a,- 



Imz\ 2 



and 

(3.2.19) £; 2 = {yeR 2rl :Ky) = 0}. 

The next lemma gives a first inclusion between these two sets E\ and E 2 . 

Lemma 3.2.1. We have 

(3.2.20) E x C E 2 . 

Proof of Lemma 3.2.1. Let y G If y = then y belongs to E 2 since from (3.2.15), 
p is a quadratic form in the variable y. If y ^ 0, we set 

2n 



E2 iv/ /Re z 

> and Vj = 1, 2n, % = y — — 



We recall from (3.2.12) that z e B\ d'E(q) implies that Re z > 0. Then, since, on 
one hand 

2n 

^y 2 = Rez, 

3 = 1 
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and that, on the other hand, we have from (3.2.17) and (3.2.18) that 

2n -~ 2n . . 2n T 

E_ 2 Ke z v-^ 2 tie z ^— % 1m 2 2 

j=i i =1 i =1 

because y <S £a, we deduce from (3.2.16) and the homogeneity of degree 2 of the 
quadratic form p that 

,_. Re z 

p(y) = —j—p(y) = o. 

According to (3.2.19), this proves that y E E 2 and ends the proof of the lemma 3.2.1. □ 

Then, we can notice from (3.2.14) that the boundary of the numerical range <9£(g) 
is given by 

(3.2.21) (1 + ia 1 )R + U (1 + ia n )R + . 
Since the numerical range X(q) is a closed set, the assumption 

zeB\9%)cE(g)\aE( (( ) = %), 
induces from (3.2.21) that 

Im z 

This implies that the signature (ri,Si) of the quadratic form r defined in (3.2.18) 
fulfills 

(3.2.22) (n, si) efxf and n + s 1 < 2n. 
Thus, we can assume after a new labeling that 

(3.2.23) r(y) = a±y\ + ... + a ri y^ - a ri+1 yl i+l - ... - a ri+si y* i+si , 

with a,j > for all j = 1, r\ + Si. It follows from (3.2.17) and (3.2.23) that in these 
new coordinates, the set E\ is the direct product of a proper cone C of R ri+Sl and 

(3.2.24) Ei=CxR 2 "- ri - Sl . 



Figure 8 




xR 2n-ri- S i _^ R 2n^ 



c 

We are now going to prove that the two sets E\ and E 2 are equal 

(3.2.25) E l =E 2 . 

Let us reason by the absurd by assuming that it is not the case. Then, we could find 
from the lemma 3.2.1, 

(3.2.26) y eE 2 \E 1 ,y = (y' ,y'^ with y' E W^+ s \ $ E R^- ri -s K 

We deduce from (3.2.24) that y' £ C. Let us now recall an elementary geometrical 
fact that we will use several times. This fact is that the intersection of a real line and 
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a real quadric surface is reduced to either 0, 1 or 2 points, or the line is completely 
contained in the quadric surface. We first begin by proving that 

(3.2.27) R ri+Sl x {y" = y'^ C E 2 . 

Indeed, let us consider the affine subspace 

F = {y e K 2 " : y = (y',y") G K ri+Sl X R 2n ~^~ s \ y" = y»}. 

We identify for more simplicity the space F to the space M ri+Sl . We agree to say that 
a point x' of K ri+Sl belongs to the set E 2 to mean that the point (x' ,y'o) belongs to 
the set E 2 . With this convention, it is sufficient for proving the inclusion (3.2.27) to 
consider some particular lines of R ri+Sl , containing the point y' defined in (3.2.26) 
and, which have an intersection with the cone C in at least two other different points 
u' and v' (see Figure 9). These lines are necessary contained in the quadric surface 
E 2 because from the lemma 3.2.1, 

Ei C E 2 , 

and that there are at least three different points of intersection between these lines 
and the quadric surface E 2 , 

(«0. V'o) e C X =EiC E ^ (u / ; ^ ^n-rr-s, = ^ c ^ 

and (y' , y'o) & E 2 . Thus, we prove that the shaded disc appearing on the figure 10 
is completely contained in the set E 2 . By using the cone structure of the set E 2 , 
we can deduce that all the interior of the cone C (see Figure 11) is contained in E 2 . 
Then, using again other particular intersections with some lines as on the figure 12, 
we deduce from our identification of the space F to ]R ri+Sl that the inclusion (3.2.27) 
is fulfilled. 

Figure 9 



We now prove that under these conditions, we have the identity 
(3.2.28) E 2 =M. 2n . 

Indeed, let us consider (#,,$,') G K 2 " = M ri+S1 x l^n-n-^ If ^ g ^ thcn 

(y'o,y'o) e£ 2 , 
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Figure 10 

C 




C 



because from (3.2.20) and (3.2.24), (y' ,y'o) € E 1 and E 1 C E 2 . If, on the other hand 
y' £ C, we can choose a point u e M ri+;Sl different from y' such that u £ C, and such 
that the line containing y' and u in M ri+Sl , has an intersection with C in at least two 
other different points v and w (see Figure 13). Thus, we can find some distinct real 
numbers t\, t<i £ M \ {0, 1} such that 

v = (1 - ti)j/o + iiu e C and w = (1 - t 2 )y' + t 2 u G C. 
Considering now the line 

D = {(l-t)(jj' ,yZ)+t{u,yZ):teR}, 
we can notice that this real line contains at least three different points of E 2 : 

(v^l-t^+hy'J), (w,(l-t2)y%+t 2 yZ) and (u,i#)- 

Indeed, this is a consequence of the fact that v and w belong to C, and from (3.2.20), 
(3.2.24) and (3.2.27). Thus, the line D is contained in the quadric surface E 2 . This 
implies that {y' ,yo) efc E 2 . 

To sum up, we have proved that if the two sets E\ and E 2 are different then the 
set E 2 is equal to M 2 ". This fact induces in view of (3.2.19) that the quadratic form p 
is identically equal to zero. By coming back to the first coordinates (x, £) = Py, it 
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Figure 12 




C 



follows from (3.2.15) that the quadratic form {Re g,Im q} is also identically equal to 
zero, which contradicts (3.2.2). This proves the identity (3.2.25), 

Ei = £/2- 

With this fact, we can resume our first reasoning by the absurd, which assume in 
(3.2.11) the existence of a point z e B \ dY,{g). Let us now consider yo E\ = E^. 
This is possible according to (3.2.2), (3.2.15) and (3.2.19). We deduce from (3.2.17) 
and (3.2.19) that r(y ) and p(yo) are non-zero. By considering AgK* such that 

p(yo) = Ar(yo) 

and 



(3.2.29) 



r(y) =p{y) - Ar(y), 
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it follows from (3.2.17), (3.2.19), (3.2.25) and (3.2.29) that 

(3.2.30) E 1 C {y G K 2 " : f(y) = 0}. 
This inclusion (3.2.30) is strict since 

f(y ) = and y E x . 

By using now exactly the same reasoning as the one previously described to prove 
(3.2.25), about the intersections of real lines and quadric surfaces, we prove that the 
quadratic form f is necessary identically equal to zero. Then, it follows from (3.2.29) 
that 

(3.2.31) p = Xr. 

By coming back to the first coordinates (cc, £) = Py, we get using (3.2.14), (3.2.15), 
(3.2.18) and (3.2.31) that for all (x,£) G R 2, \ 

(3.2.32) {Re q, Im q}(x, £) = \ (im q(x, £) - — Re q(x, Q) . 

Let us now consider (xo,£o) G R 2 ™ such that q(xo,£o) € 9T,(q) \ {0}. This is possible 
since the numerical range X(g) is a closed angular sector with a top in and a positive 
opening. We deduce from (3.2.5) and (3.2.8) that we necessarily have 

{Re q,Im <?}(x ,£ ) = 0. 

This induces from (3.2.32) that 

(3.2.33) Imq(x ^ ) = l ^Req(x a ,Z ), 

Re z 

because A G K*. Since according to the shape of the numerical range and 
(3.2.12), 

<z(x ,£ ) G 0£(«) \ {0} C {z G C : Re z > 0}, 

the identity (3.2.33) proves that the point z also belongs to the set dT,(q), but it 
contradicts the initial assumption 

z G B\dH(q). 

Finally, this ends our reasoning by the absurd and proves (3.2.6). 

3.2.1.b. Existence of semiclassical quasimodes at the interior of the numerical range. 
To prove the existence of semiclassical quasimodes for the associated semiclassical 
operator 

(q(x,hO W )v<h<l, 

in every point of the numerical range's interior (Theorem 2.2.1), we use an existence 
result of semiclassical quasimodes for general pseudodifferential operators violating 
the condition ('I') 5 . Let us mention that this result generalizes the two existence 
results of semiclassical quasimodes given by E.B. Davies, in the case of Schrodinger 
operators (Theorem 1 in [4]), and by M. Zworski in [17] and [18], for pseudodifferential 
operators. 

This existence result of semiclassical quasimodes can be stated as follows. Let us 
consider a semiclassical symbol P(x, £; h) in S(((x, £)) m , dx 2 + c?£ 2 ) with m G R+, 

((x,C)) 2 = l + x 2 +e 2 , 



'The definition of the condition ('I') is recalled below. 
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where S (((x , ^)} m , dx 2 + d£ 2 ) stands for the following symbol class 
S(((x,(,)) m ,dx 2 +d£ 2 ) = {a(x,£;h) e C°°(R£ x R r £ , C) : 

Va G N 2 ", sup ||<(x,0)- m a?, c a(a:^;/»)|Uco (Ra » ) <+ooj, 
with a semiclassical expansion 

+ OC 

(3.2.34) P(x,£;h)~J2 hj P^ x >®' 

3=0 

where for all j E N, pj is a symbol of the class S(((x, £)) m , dx 2 + d£, 2 ) independent 
from the semiclassical parameter h. 

Let z E C, we assume that there exists a function q € C£°(R 2 ™, C), where 
C£°(R 2 ™,C) stands for the set of bounded complex- valued functions on R 2 ™ with 
all derivatives bounded, and a bicharacteristic curve, t E [a, 6] i— > 7(i), of the real part 
Rc(<7o(po — z )) OI the symbol qo{po ~ z), with a < 6, such that 

(3.2.35) Vt e [a,b], q a ("f(t)) ^ and 

Im [q (7(a) )(po (7(a)) - «)] > > Im[ go (7(&))(Po(7(&)) - *)] • 

Theorem 3.2.1. Under these assumptions (3.2.34) and (3.2.35), /or aiZ open neigh- 
bourhood V of the compact set 7([a, b]) in R 2 ™ and for all N E N, there exist ho > 
and (uh)o<h<h a semiclassical family in S(R n ) such that 

\\uh\\L2(Rn) = 1, FS((u h ) 0<h <h„) C V and \\P(x, h£; h) w u h - zu h \\ L 2 {Rn) = 0(h N ), 
when h — > + . 

The notation FS((u/ l )o</ l </ lo ) stands for the frequency set of the semiclassical fam- 
ily {uh)o<h<h defined as the complement in R 2 ™ of the set composed by the points 
{xq, Co) G R 2 ", for which there exists a symbol Xo(x, £; h) E 5(1, dx 2 + o!£ 2 ) such that 

Xo(zo,£o;/i) = 1 and ||xo(^, M^^IU^R") = 0(/i°°), 

when ft, — > + . 

This existence result of semiclassical quasimodes is an adaptation in a semiclassical 
setting of the proof given by L. Hormander in [7] for proving that the condition ($) is a 
necessary condition for the solvability of a pseudodifferential operator (Theorem 26.4.7 
in [7]). The existence of this result has been first mentioned in [5]. A complete proof of 
this adaptation in a semiclassical setting is given in [11]. This result shows that when 
the principal symbol po — z of the symbol P — z violates the condition ('J'), there exists 
in this point z some semiclassical quasimodes inducing the presence of semiclassical 
pseudospectrum of infinite index for the semiclassical operator P(x, h) w . 

Condition (*). A complex- valued function p e C°°(R 2 ™, C) fulfills the condition (*) 
if there is no complex- valued function q E C°° (R 2n , C) such that the imaginary part 
Im(gp) of the function qp changes sign from positive values to negative ones along 
an oriented bicharacteristic of the symbol Re (op) on which the function q does not 
vanish. 

By using the characterization given in the previous section for the interior of the 
numerical range S(g) (see (3.2.4) and (3.2.6)), we are now going to prove that the 
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principal symbol q(x,£) — z of the scmiclassical operator 

q(x,hO w -z, 

violates the condition for all z in T,(q). This violation of the condition (^) will 
induce in view of the theorem 3.2.1 that for all z € £((?) and N € N, we can find a 
semiclassical quasimode {uh)o<h<h a € with /i > 0, verifying 

H^Hl^R") = 1 and \\q(x, h^) w u h - zUh\\^(M.n) = 0(h N ) when h -> + , 
which will end the proof of Theorem 2.2.1. 

Let us consider z € £(<?). We are now going to prove that there is actually a 
violation of the condition (<&) for the symbol q — z. According to (3.2.4) and (3.2.6), 
there are two cases to separate. 

Case 1. Let us assume that there exists (xo,Co) G such that 

(3.2.36) z = q(x , Co), {Rc(<7 - z), lm(q - z)}(x , Co) = {Re q, Im <?}(:z;o, Co) < 0. 
By considering the solution of the following Cauchy problem 

we define the following function 

(3.2.38) f(t) = Im q(Y(t)) - Im g(x , Co)- 

As mentioned before, (3.2.37) is a linear Cauchy problem. It follows that its solution 
Y is global and that the function / is well-defined on R. A direct computation using 

(3.2.37) and (3.2.38) gives that for all t G R, 

(3.2.39) f(t) = {Req,lmq}(Y(t)). 
Since from (3.2.36), (3.2.37), (3.2.38) and (3.2.39), 

/(0) = 0, /'(0) = {Re g,Im g}(x ,Co) < 

and i?R e 9 -Ro 2 = Hr c q , we deduce in this first case that the imaginary part of the 
function q — z changes sign, at the first order, from positive values to negative ones 
along the oriented bicharacteristic Y of the symbol Re q — Re z. This proves that the 
symbol q — z actually violates the condition (\I>). 

Case 2. Let us now assume that there exists (xo,£o) € K 2 ™ such that 

(3.2.40) z = q(x ,£ ), {Re(q - z), Im(q - z)}(x , Co) = {Rc q, Im q}(x , Co) > 0. 

We consider as in the previous case, the global solution Y of the Cauchy problem 
(3.2.37) and the function / defined in (3.2.38). Since from (3.2.37), (3.2.38), (3.2.39) 
and (3.2.40), 

(3.2.41) /(0) - 0, /'(0) = {Re q, Im q}(x , Co) > 0, 

we deduce this time that the imaginary part of the function q — z also changes sign, 
at the first order, along the oriented bicharacteristic Y of the symbol Re q — Re z. 
Nevertheless, this change of sign is done in the "wrong" way. It is a change of sign 
from negative values to positive ones, which does not induce directly a violation of 
the condition (<&). To check that there is actually a violation of the condition ("J") 
in this second case, we need to study more precisely the behaviour of the function 
Im q — Im z along this bicharacteristic Y. 

We deduce from (3.2.41) that there exists e > such that 

vte [-e.e], f(t) >o, 
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which induces that 

(3.2.42) /(e) > and /(-«•) < 0, 

since from (3.2.41), /(0) = 0. By using the following lemma, we obtain that for all 
S > 0, there exists a time to(S) > e such that 

(3.2.43) \Y(t (S))-Y(-s)\<5. 



Figure 14 




Lemma 3.2.2. IfY(t) = (x(t),£(t)) is the C°°(R,K 2 ™) function solving the linear 
system of ordinary differential equations 

Y'(t) = H neq {Y(t)), 

where Re q is the symbol defined in (3.2.7), then we have 

Vt G R, Ve > 0, VM > 0, 37\ > M, 3T 2 > M, 

\Y(t ) - Y(t + Ti)| < e and \Y(t ) - Y(t - T 2 )| < e. 



Proof of Lemma 3.2.2. If Y(t ) = (a 1} a n , b\, b n ) G R 2 ™, we deduce from (3.2.7) 
that the function Y(t) = (x (£),£(£)) solves the following Cauchy problem 

f x' j (t) = 2\ j £ j (t) 
Cj(t) = -2X jXj (t) 



Vj = l,...,n, { 



Xj(t ) 

{ &(*o) = 



6,, 



(3.2.44) 



It follows that for all j — 1, n and t G R, 

£j(t) = 6j sin (2(i — io)Aj) + aj cos (2(i — io)Aj) 
£j(t) = bj cos (2(i - t )\j) - Oj sin (2(t - i )Aj) . 

Setting f3j = Xj/tt for all j = 1, ...,n, we need to study two different cases. 

Case 1: Vj G {l,...,n}, /3j G Q. In this case, the function Y is periodic and the 
result of Lemma 3.2.2 is obvious. 
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Case 2: {(3\, (3 n ) £ Q™. In this second case, we use the following classical result of 
rational approximation: Ve > 0, V(0i, 9„) G R n \ Q", 3pi, ...,p n g Z, 3q g N* such 
that 

Pj_ 
Q 



e 

< -. 

q 



< sup 

If < ei < 1/2, we can therefore find some integers ...,pi.n G Z and <? ei g N* 
such that 

0< sup \q El 0j-pi,j\<e!. 

j=l,...,n 

If 

£2 = x sup |g £l /3j — pi,j| > 0, 

J — l,...,n 

using again this result of rational approximation, we can find some other integers 
P2.1, ■■■,P2,n g Z and g £2 g N* such that 

< sup \q E2 pj -p 2 ,j\ < e 2 - 

j=l,...,n 

By using this process, we build some sequences (p m ,j)meN* of Z for j = l,...,n, 
(£ ra ) m£ M* of M^J_ and (g em ) m6 N* of N* such that for all m > 2, 

(3.2.45) 0< sup \q £m Pj - p m ,j\ < e m = ]- sup \q Em _ 1 f3 J - p m - h] \ 

j—l,...,n a j — l,...,n 

and 

(3.2.46) < e rn < ^ £l . 

The elements of the sequence (q em ) m£ ^* are necessary two by two different. Indeed, 
if q ek = q £l for k < I, this would imply according to (3.2.45) and (3.2.46) that 

Vj = l,...,n, \p kJ -pi,j\ < \q ek (3j ~Pk,j\ + \q Sl Pj ~Pl,j\ <s k + ei<l, 

because < e\ < 1/2, which would induce that Vj = 1, ...,n, p k .j = pij because p k .j 
and pi.j are some integers; and would contradict (3.2.45) because 

< sup \q Sl pj -pij\ <si< - sup \q £k f3j - p k ,j\- 

j=l,...,n ^ j — l....,n 

Since the sequence (<7 £m )meN* is composed of integers two by two different, we can 
assume after a possible extraction that q €m — ► +oo when m — > +oo. We deduce from 
(3.2.44), (3.2.45) and (3.2.46) that 

Y (to + <7e m ) -» ^(*o) when m —> +oo. 

Then, considering (/3i, /3„) = (— (3i, —(3 n ), we obtain by using the same method 
a sequence (<7 £m )meN* of integers such that q £m — > +oo and 

F(to — <?e m ) — » ^(*o) when m — > +oo. 

This ends the proof of Lemma 3.2.2. □ 

Since from (3.2.42), /(-e) < 0, we deduce from (3.2.38) and (3.2.43) that there 
exists to > e such that /(to) is arbitrarily close to f(—e). It follows in particular that 
we can find t > e such that f(to) < 0. Since from (3.2.42), /(e) > and f(t ) < 0, 
we deduce from (3.2.38) and (3.2.40) that the function 

1 1 ► Im q(Y(t)) — Im z, 

changes sign from positive values to negative ones on the interval [e, to]- This proves 
that the imaginary part of the function q—z actually changes sign from positive values 
to negative ones along the oriented bicharacteristic Y of the symbol Re q — Re z; and 
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that the symbol q — z also violates in this second case the condition (W). This ends 
the proof of Theorem 2.2.1. 

3. 2. I.e. Another proof for the existence of semiclassical quasimodes. In the following 
lines, we give another proof for the existence of semiclassical quasimodes in some 
points of the numerical range's interior. The result proved in this section is weaker 
than the one given by the theorem 2.2.1, since we prove the existence of semiclassical 
quasimodes in every point of the numerical range's interior without a finite number 
of particular half-lines. 

Let us consider a non-normal elliptic quadratic differential operator 

(3.2.47) q(x,O w : B -> L 2 (R n ), 

in dimension n > 2. We assume, as before, that (3.2.7) is fulfilled. Using that 
the quadratic form Re q is positive definite, we can simultaneously reduce the two 
quadratic forms Re q and Im q by choosing an isomorphism P of R 2 ™ such that in 
the new coordinates y — P^ 1 (x,( i ), 

In 2n 

(3.2.48) r 1 {y)=Req{Py) = Y^yl r 2 (y) = Im q{Py) = ^a^ 2 , 

with ai < ... < a n . Let us study when the differential forms dri(y) and dr 2 {y) are 
linearly dependent on R i.e. when there exist (A, [i) e R 2 \ {(0, 0)} such that 

(3.2.49) Xdn(y) + ^dr 2 (y) = 0. 

It follows from (3.2.48) and (3.2.49) that for all j = 1, ...,2n, 

(3.2.50) (X + fiaj)yj = 0. 

If y ^ 0, then there exists jo <S {1, 2n} such that j/ JO ^ 0. This implies that 

(3.2.51) \ + na jo =Q. 

We deduce from (3.2.50) and (3.2.51) that yj = if ctj ^ ctj . Thus, we obtain that if 
zet(g)\((l + Ml )i;u...U(l + w n )i;), 

then the differential forms dRc q and dlvaq arc linearly independent on R in every 
point of the set q^ 1 (z). 

Figure 15 

4 

(1 + ia n )W + 
(1 + ia n ^)Rl 



(1 + ia 3 )R* + 



(1 + ia 2 )R* + 
(1 + ia 1 )R\. 
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Let us consider such a point 

z e t(q) \ ((1 + i ai )R* + U ... U (1 + ia n )R* + ) . 

Since the dimension n > 2, we can apply the lemma 3.1 in [5] (see also the lemma 8.1 
in [9]). It follows that for any compact, connected component T of q~ 1 {z), we have 

(3.2.52) ^{Rc q,Im q}{p)\ q . z {dp) = 0, 

where A g , z stands for the Liouville measure on q~ 1 (z), 

A„ z A dRe q A dim q = — - . 

n\ 

The set q~ x (z) is a non-empty submanifold of codimension 2 in R 2 ™. We deduce from 
(3.2.4) and (3.2.6) that there exist (xq,£o) <= 9 _1 ( z ) sucri that 

(3.2.53) {Re g ,I m(Z }(xo,eo) ^ 0. 

Then, it follows from (3.2.52) and (3.2.53) that there necessary exists (xq, £o) G <7 _1 (-2) 
such that 

(3.2.54) {Re g,Im «}(£,,, fo) < 0. 

Under this condition (3.2.54), we can use the reasoning given in the first studied case 
(see (3.2.36)) to prove that the imaginary part of the function q—z changes sign, at the 
first order, from positive values to negative ones along an oriented bicharacteristic of 
the symbol Re q — Re z. This induces that the symbol q—z violates the condition (\I>); 
and we can conclude by using the theorem 3.2.1. Let us mention that we can also 
directly use the existence result of semiclassical quasimodes given by M. Zworski in 
[17] and [18]. This second proof gives the existence of semiclassical quasimodes in 
every point belonging to the set 

±(q) \ ((1 + ian)R* + U ... U (1 + ia n )M* + ). 

3.2.2. On the pseudo spectrum at the boundary of the numerical range. In this section, 
we give a proof of the theorem 2.2.2. Let us consider a non-normal elliptic quadratic 
differential operator 

q(x,O w :B^L 2 (R n ), 

in dimension n > 1. We assume that S(g) ^ C, and that its Weyl symbol q(x, £) is of 
finite order kj on a half-line Aj, j € {1, 2} (See the definition given in (2.2.9)), which 
composes the boundary of its numerical range 

(3.2.55) dY,(q) = {0} U A x U A 2 . 

As we have already done several times, we can reduce our study to case where (3.2.7) 
is fulfilled. 



Proof of Theorem 2.2.2. Let us consider the following symbol belonging to the 
C£°(R 2 ™,C) space, composed of bounded complex- valued functions on R 2n with all 
derivatives bounded 

(3.2.56) r(*,0- q ^- Z 



i + x 2 + e ' 



with z G Aj. Setting S(r) = r(R 2 ™), we can first notice that 

z e 8T,(q) \{0}^0e dt(r). 
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Let us also notice that the symbol r fulfills the principal-type condition in 0. Indeed, 
if (xo,^o) € R 2 ™ was such that r(xo,Co) = and dr (xo,t;o) = 0, we would get from 

(3.2.56) that 

(3.2.57) dq(x ,£ ) = 0. 
Since from (3.2.7) and (3.2.57), we have 

n 

dRe q(x ,S,o) = 2 Aj {{x )jdxj + {£o)jd£j) = °> 
j=i 

this would imply that 

(xo,Zo) = (0,0), q(x ,Zo) = 0, 

because q is a quadratic form and that Xj > for all j = 1, n. On the other hand, 
since r(xo,^o) = 0, we get from (3.2.56) that g(xo,£o) = z^0 because 

z e A 3 c as(«) \ {0}, 

which induces a contradiction. It follows that the symbol r actually fulfills the 
principal-type condition in 0. Let us notice that, since symbol q is of finite order kj 
in z, this induces in view of (3.2.56) that the symbol r is also of finite order kj in 0. 
On the other hand, we deduce from (3.2.7) and (3.2.56) that the set 

{(x,Q G R 2 " : r(x,0 = 0} = {(x,£) G M 2 " : q(x,Z) = z}, 

is compact. Under these conditions, we can apply the theorem 1.4 in [5], which proves 
that the integer kj is even and gives the existence of positive constants ho and C\ 
such that 

(3.2.58) V0<fe</i o ,V«€ S(R n ), \\r(x, h£) w u\\ L 2 {Rn) > dh^ +1 \\u\\ L 2 (M _ny 



Remark. We did not check the dynamical condition (1.7) in [5], because this assump- 
tion is not necessary for the proof of Theorem 1.4. Indeed, this proof only use a part 
of the proof of lemma 4.1 in [5] (a part of the second paragraph), where this condition 
(1.7) is not needed. 



By using some results of symbolic calculus given by Theorem 18.5.4 in [7] and (3.2.56), 
we can write 

(3.2.59) r(x, h£) w {l + x 2 + h 2 (, 2 ) w = q(x, h£) w -z + hn{x, h£) w + h 2 r 2 (x, h£) w , 
with 

(3.2.60) ri (x, = -ix^(x, + i^(x, 
and 

(3-2.61) r 2 (x,0 = - — (*,€) 

We can easily check from (3.2.56) that these functions r\ and r 2 belong to the space 
C£°(R 2 ™, C), and we deduce from the Calderon-Vaillancourt theorem that there exists 
a positive constant C 2 such that for all u E S(M. n ) and < h < 1, 

(3.2.62) \\n(x,hO w u\\ L 2 < C 2 \\u\\ L 2 and \\r 2 {x,K) w u\\ L 2 < C 2 \\u\\ L 2. 
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It follows from (3.2.58), (3.2.59), (3.2.62) and the triangular inequality that for all 
u e 5(R") and < h< ho, 

Ci^||(l + a; 2 + ftVr«l|L"(R») 

< \\r{x,hira+x 2 +h 2 eru\\ LHm 

< \\q(x, h£) w u - zu\\ L 2( Rn ) + C 2 h(l + h)\\u\\ L 2 {Rn) . 

Since from the Cauchy-Schwarz inequality, we have for all u G S(M. n ) and < h < 1, 
||«Hi2(R») < ll«lli2(R») + \\xu\\ 2 L 2 {Kn) + ||/li3xU||i2 (R „ ) 

= ((i+x 2 +h 2 e) w u,u) L2{Rn) 
< \\(i+x 2 +h 2 e) w u\\ L 2 {Rn) \\u\\ L 2 {Rn) , 

we obtain that for all u e 5(M") and < h < ho, 

(3.2.63) Cift^||«||£2(R») < ||g(a;,^) t, '«-^llL2( R ")+C 2 ^(l + /i)|| U || i 2 (R „ ) . 

Since fcj > 1, we deduce from (3.2.63) that there exist some positive constants h' and 
C 3 such that for all < h< h' and ue5(M n ), 

||9(^ 5 - OT lli 2 (l") > C 3 /l^ + 1 ||u||i2( M n). 

Using that the Schwartz space 5(K") is dense in B and that the operator 

q{xM) w +z, 

is a Fredholm operator of index 0, we obtain that for all < h < h' , 

\\(q(x,hO w -z)~ 1 \\<C3 1 h-& 1 , 
which ends the proof of Theorem 2.2.2. □ 

About the case of infinite order, the situation is much more complicated. As 
mentioned before, we cannot expect to prove a stronger result than an absence of 
semiclassical pseudospectrum of index 1, but we can actually prove that there is never 
some semiclassical pseudospectrum of index 1 on every half-line of infinite order, by 
using a result of exponential decay in time for the norm of contraction semigroups 
generated by elliptic quadratic differential operators proved in [12]. 

The result proved in [12] shows that the norm of a contraction semigroup 

lle'^Hk^), t>0, 

generated by an elliptic quadratic differential operator q(x,l;) w with a Weyl symbol 
verifying 

Reg<0, 3(^0,^) eK 2 ™, Re q{x ,£o) ± 0, 
decreases exponentially in time 

(3.2.64) 3M,a > 0,Vt > 0, \\ e tq(x ^ \\ C(L 2 } < Me~ at . 

Let us consider a non-normal elliptic quadratic differential operator 

q(x,O w :B^L 2 (R n ), 

in dimension n > 1 such that ^ C. We explain in the following lines how (3.2.64) 
allows to prove that there is never some semiclassical pseudospectrum of index 1 on 
any open half-lines composing the boundary of the numerical range dT,(q) \ {0}. 
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Let z e dT,(q) \ {0}. Since the numerical range is a closed angular sector with 
a top in and a positive opening strictly lower than n, we can find e € {±1} such 
that 

(3.2.65) Rcieiz^q) < 0, 3(x ,£o) € K 2 ", Re(£i2 _1 g)(a;o, £o) 7^ 0. 
Using the theorem 2.8 in [2], we obtain that for all rj g M, 

= - w~ 1 e(ei»?-e»«~ 1 g(x,0 , ")~ 1 

(3.2.66) = -irt / e^V 6 ** '^^"ds. 

Jo 

It follows from (3.2.64) and (3.2.65) that for all r\ e M, 

r+OO 

IKg^O' -^) || < M" 1 / l|e s£i2 " 1<?(x ' ?r |U(^)^ 

Jo 

< Izl" 1 / Me~ as ds = Ur 1 — < +oo, 
Jo a 

which proves the absence of semiclassical pseudospectrum of index 1 on the half-line 

zR+. We can actually use the theorem 2.8 in [2] because 

tKc C\a(eiz- 1 q(x,£,) w ). 

Indeed, if it was not the case, we would deduce from (2.1.7) that there exists u € 
B \ {0} and A G R such that 

eiz^ 1 q(x, S,) w u = i\o u o- 

Since from (3.2.65), the quadratic form — Re(eiz~ 1 q) is non-negative, we deduce from 
the symplectic invariance of the Weyl quantization and the theorem 21.5.3 in [7] that 
there exists a metaplectic operator U such that 

k k+l 

(3.2.67) -[Re{eiz-\{x,0)] W = U-'(Y j \ j {Dl j +x>)+ ]T x])u, 

j=l j=k+l 

with k, I € N and Xj > for all j = 1, k. By using that U is a unitary operator on 
L 2 (R n ), we obtain that 

= - Rc(i\ u ,u ) L 2 

= -Re(eiz~ 1 q(x,^) w u ,u ) L2 

= - ([Re(eiz~ 1 q{x,£))] w u ,u ) L2 

k k+l 

= E x i(\\ D *i Uu o\\h + n^«oiii0+ E \\ x i Uu °\\h, 

j = \ j=k+l 

which induces that uq = 0, because from (3.2.65) and (3.2.67), k + I > 1. It follows 
from (2.1.7) that there exists £o > such that 

<j(eiz- 1 q(x,£) w ) C {z e C : Re z < -e }. 
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